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Abstract 

We investigate quantum entanglement of a scalar field in the inflationary universe. By introduc- 
ing a bipartite system using a lattice model of scalar field, we apply the separability criterion based 
on the partial transpose operation and numerically calculate the bipartite entanglement between 
separate spatial regions. We find that the initial entangled state becomes separable or disentangled 
when the size of the spatial regions exceed the Hubble horizon. This is a necessary condition for 
the appearance of classicality of the quantum fluctuation. We further investigate the condition for 
the appearance of the classical distribution function and find that the condition is given by the 
inequality for the symplectic eigenvalue of the covariance matrix of the scalar field. 
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I. INTRODUCTION 



Inflation provides a mechanism of generation of primordial inhomogeneity which is needed 
for the formation of large scale structures in our present universe. During the qusi-de Sitter 
expansion stage of the inflationary universe, short wavelength quantum fluctuations of the 
inflaton field are generated by particle creations and then they are stretched by the cos- 
mic expansion and their wavelength exceeds the Hubble horizon beyond which the physical 
process proceeds independently. After the wavelength of generated quantum fluctuations 
become larger than the Hubble horizon, the quantum nature of the fluctuation is expected 
to be lost and the statistical property of fluctuations is replaced by the classical distribu- 
tion function. This is the assumption of the quantum to classical transition of the quantum 
fluctuation generated during the inflation. Once this assumption is adopted, we can use gen- 
erated "classical" fluctuations as initial perturbations for the large-scale structure formation, 
which obeys deterministic classical dynamics. 

We must explain or justify this assumption of the quantum to classical transition of 
riraordial fluctuation and many investigations have been done on this subject so far|l|, y, 

111 . 112] . A rough outline of the classicalization mechanism proposed 
M\ is as follows: In the inflationary universe, particle creations 
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occur due to the accelerated expansion of the Universe and the quantum field becomes a 
highly squeezed state. For such a highly excited state as a squeezed state, noncommutativity 
between canonically conjugate variables can be neglected and the operator nature of variables 
are effectively lost when we evaluate the expectation value of quantum variables. This means 
we can regard operators as c-numbers. Furthermore, it can be shown that there appears a 
sharp peak around a line of the Wigner function for the state of quantum fluctuations and 
this indicates establishment of classical correlation between canonically conjugate variables 
on the phase space. At this stage, the Wigner function itself can be interpreted as a classical 
probability distribution function and the quantum fluctuations themselves can be treated 
as classical stochastic variables. Hence, after sufficient squeezing, we can represent the 
nature of the quantum fluctuations by the classical stochastic variables with an appropriate 
probability distribution function of which property is determined by the state of the quantum 
fluctuations. For the classicality of the quantum system, we also need the mechanism of 
decoherence and the studies along this line have been done with the assumption of an 
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appropriate coupling between the system and the environment [2], y, 

However, these analysis do not pay attention to an important aspect of quantum me- 
chanics, quantum entanglement, which definitely distinguishes a quantum world from a 
classical one. When we calculate a correlation function of observables between spatially 
separated two regions, we have a possibility that the quantum correlation function cannot 
be reproduced using a local classical probability distribution function if these two regions 



are entangled 



13 



14| and the classical locality is violated. In other words, we cannot regard 
the quantum fiuctuations as the classical stochastic fluctuations as long as the system is 
entangled. Therefore, it is important to clarify the relation between the entanglement and 
the appearance of the classical nature to fully understand the mechanism of the quantum 
to classical transition of primordial fluctuations. 

In this paper, we consider the entanglement of quantum fleld between two spatially sepa- 
rated regions in a de Sitter universe and aim to understand the mechanism of the quantum to 
classical transition of fluctuation from the viewpoint of quantum entanglement. This paper 
is organized as follows: We introduce the concept of entanglement and separability in Sec. II. 
In Sec. Ill, we calculate the bipartite entanglement of a scalar fleld in the de Sitter universe. 
In Sec. IV, we discuss the relation between the entanglement and the quantum to classical 
transition in the inflationary universe. Sec. V is devoted to summary and conclusion. We 
use units in which c = h = SnG = 1 throughout the paper. 

II. SEPARABILITY AND ENTANGLEMENT 

In this paper, we focus on a bipartite system composed of two Gaussian modes described 
by canonical variables {qi,Pi) and {q2,P2)- This is the simplest case deriving entanglement 
for continuous variable system. Let Alice be in possession of mode 1 and let Bob be in 
possession of mode 2. A quantum state p of the bipartite system is deflned to be separable 



if and only if p can be expressed in the form 15 1 



j j 

where pjA and pjB are density operators of the modes of Alice and Bob, respectively. If 
the state of the system cannot be expressed in this form, the quantum state of the system 
is called entangled. When the state is entangled, the observables associated to the party 
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A and B are correlated and their correlations cannot be reproduced with purely classical 



means. This 



eads to the phenomena peculiar to th e q uantum mechanics such as the EPR 



correlation jl3|] and the violation of Bell's inequality 14l |. 



For a bipartite Gaussian two mode system, we have necessary and sufficient conditions 
for separability and we can judge whether the system is entangled or not using these 



criteria. In this paper, we adopt a criterion proposed by Simon 16^ which uses the partial 
transpose operation for a bipartite system. We define the phase space variables as 

Pi 

Using these variables, the canonical commutation relations are expressed as 



(2) 



a,/5= 1,2,3,4, 
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J 




(3) 



The Gaussian state is completely characterized by the covariance matrix 



Vap = ^ (^iaip + ida) = ^tr ( {^0.^13 + ^p^a)p 



(4) 



where we assume the state with (^q,) = 0. For a physical state, the density matrix must be 
non-negative and the corresponding covariance 



V + -O > 

2 - 



(5) 



which is the generalization of the uncertainty relation between two canonically conjugate 
variables. The separability of the bipartite Gaussian state is expressed in terms of the partial 
transpose operation defined by 



I =Ai, A = diag(l, 1,1,-1) 



(6) 



This operation reverses the sign of Bob's momentum. With this operation, the covariance 
matrix transforms as 

V = AVA^. (7) 



The necessary and sufficient condition of the separabihty is given by the inequahty 

v + ^n>o (8) 

which represents the physical condition for the partially transposed state. For an entangled 
state, this inequality is violated and the partially transposed state becomes unphysical. To 
measure the degrees of entanglement, we introduce the logarithmic negativity via symplectic 
eigenvalues of the covariance matrix. The covariance matrix can be diagonalized by an 



appropriate symplectic transformation S G Sp(4, R), SflS = 17 as follows 



HQ 



SVS' = diag(i/+, Z/+, z/_, z/+ > z/_ > (9) 

where are symplectic eigenvalues. In terms of symplectic eigenvalues, the physical con- 
dition (j5]) can be expressed as 

^- > I (10) 
and the separability condition ([8]) can be expressed as 

i>_ > -. 11) 

The logarithmic negativity is defined by 

= -min[log2(2z>_),0]. (12) 

For an entangled state, z/_ < 1/2 and we have En > 0. For a separable state, z/_ > 1/2 and 
we have Ejy = 0. Practically, the symplectic eigenvalues can be obtained as eigenvalues of 
the matrix zfivHQ. 

III. ENTANGLEMENT OF QUANTUM FIELD IN THE DE SITTER UNIVERSE 
A. One-dimensional lattice model of scalar field 

To comprehend the behavior of the entanglement of quantum fields in the inflationary 
universe, we consider a real massless scalar field in the de Sitter universe. The metric and 
the Lagrangian are 

ds'^ = aiviYi—drf + dx^), a = —-pz—, —oo < t] < 0, 

Hrj 

L = [ d^x^ , (13) 



2^ 
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where t] is the conformal time and H is the Hubble parameter, which is assumed to be 
constant in time. By introducing a conformally rescaled variable q = a(f), 



, a 
q q 



and the equations of motion of the scalar field is 



—q - djq = 0, 



(14) 



(15) 



where ' denotes the derivative with respect to the conformal time rj. 

To investigate the property of the quantum entanglement of the scalar field, we adopt 
a discrete lattice model of the scalar field in this paper. This model introduces a cutoff 
of short wavelength mode of the scalar field, which regularizes the ultraviolet divergence 
of the vacuum fluctuations. The same model is used to investigate the spatial structure of 



entanglement in a Minkowski spacetime 



201 ]. To simplify the analysis, we assume that the 



scalar field depends only on one spatial coordinate and the space is one-dimensional. Then 
the lattice version of the scalar field Lagrangian is 



L 



Ax 

~9' 2-^ 



(16) 



where qj denotes the scalar field at the j-th lattice site. Ax is a lattice spacing and is the 
total number of lattice sites. The equation of motion is 



9'' 



a" 1 

—(Ij + [2?i - a{qj+i + qj-i)] = 0, 



1,2, 



(17) 



Qo — Qn, Qn+i — Qi 



where we assume a periodic boundary condition for qj, and the parameter a 7^ 1 is introduced 
to regularize the infrared divergence which appears in the correlation function of the scalar 
field. This divergence is peculiar to one-dimensional massless scalar field. The nonunity 
value of the parameter a corresponds to adding a small mass to the scalar field 

o 2(1 -a) 



m 



(18) 



(Ax)2 

and we choose the value of a sufficiently close to unity so that the our result of calculation 
does not depends on the value of this cutoff parameter. By rescaling the time variable as 
T] — > T]Ax, the equation of motion can be written as 



a 



(19) 



The Hamiltonian is 

^ ri a' 1 

H^Yl 2P^j + 1j -^QjQj-i + -Pj1j ■ (20) 
j=i L J 

To quantize this system, we introduce the Fourier expansion of the scalar field on the lattice 
as follows 

ft = ;^E*^*"'^ k-^JLpI^"''- = (21) 

^-"^ fe=0 ^ k=0 



The equation of motion for the Fourier mode Qk is 

q'k + (^^l - 7) = 0' ^fe = 2(1 - acos^fe). (22) 

Introducing creation and annihilation operators, the quantized canonical variables are rep- 
resented as follows 

N-l 

'N 



= ^ E (/'^^^ + A*«iv-.) e-'"'. (23) 

1 

= 7^ E(-^) (^'^"^ - 9l^^N-k) (24) 

fc=0 

[^j^Pt] = ^^je, = ^kiM, [ofei, Oika] = [ofcl, QfcJ] = 0, (25) 



f^+{u;l--]h^O, g,^i(f',-^f,), hfl'-fUk-i. (26) 



As the quantum state of the scalar field, we assume the Bunch-Davis vacuum state and it 
corresponds to the following form of the mode functions 

The two point correlation functions between the canonical variables on the lattice sites are 
given by 

1 1 

9\3-t\ = ^{QjQi + Mj) = ^^\fkf cos{9k{j -i)), (28) 

fc=0 

1 1 

= 2^m+PiPj) = ]7 E l^'^l^ cos{ek{i - i)), (29) 

k=0 

1 1 ^"^ i 

f^\j-e\ = ^(^jPe + Mj) ^^Yl 2^fk9l - fk9k) cos{0k{j - £)). (30) 

k=0 

Now, we define a bipartite system using this lattice model. As we are interested in 
the correlation and the entanglement between different spatial regions, we introduce the 
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following block variables by spatially averaging the variables in given regions A and B (see 
Fig. [ID . 

QA = ^y^Qj^ PA = -^y^Pj^ = ^J' = ~^YPr (31) 

^ j&A ^ jeA ^ j(^B ^ jeB 



X 



n d n 

FIG. 1: A bipartite system in our one-dimensional lattice model. 

The region A and B contain n lattice sites and the coarse-grained field values are assigned 
to each regions. The separation between A and B is d. The commutation relations between 
these coarse-grained variables are 

[u^Pa] = [qb^Pb] =i, [qA,pB]=0 (32) 

and the set of canonical variables {qA,pA, QbiPb) constitutes a bipartite system. The covari- 
ance matrix of this bipartite system is given by the following symmetric 4x4 matrix: 

'ac\ C^h'''\, (33) 

C A) {as a2 Ics Ca/ 



i,jeA 

02 = {P\) = iPl) = ^Y1 

i,3eA 

as = ^iupA+pAU) = -Yl 

^ Til 

i,jeA 

I,. . . . , 1 
Cl = T^ilAlB + qBQA) = - > 9\i-j 



C2 = \kpApB^pBpA} = - X] ^l^-^l' 

ieAjeB 

CS = T^i^ApB + PbU) = - ^2 

i<^A,j€B 



As we do not observe the degrees of freedom of outside the regions A and B, the evolution 
of this bipartite system is nonunitary. Thus, we take into account the effect of decoherence 

8 



through our definition of bipartite system. Using these components of the covariance matrix 
V, the symplectic eigenvalues are given by 

(z/_)^ = aia2 - al + C1C2 - C3 - \a1C2 + a2Ci - 20303!, (34) 
{u- f = aia2 - 03 - C1C2 + C3 - |(aiC2 - a2Ci)^ + 4(aiC3 - a3Ci)(a2C3 - 0302)]^^^, (35) 

and we can apply the separability criterion (|TT|) to judge whether the system is separable or 
entangled. 

B. Numerical result 

We calculated the logarithmic negativity of this system numerically. The number of 
lattice sites is = 100 and the value of the infrared cutoff parameter is chosen to be 
a = 0.9999. 




FIG. 2: The dependence of the separation d on E]\[ for region size n = 4. The left panel is E'a? at 
?7 = — 10 and the right panel is at ij = —0.9. 

Fig. [2] shows the logarithmic negativity En as a function of separation d between the regions 
A and B with the region size n = 4. Initially (r/ = —10, the left panel), 7^ for d = 0, 
and En = for d > 1. The regions A and B are entangled for d = and separable for d > 1 
. This implies that the system is intrinsically entangled at this time because the choice of 
the separation d corresponds to the choice of measurement; how to observe the system. As 
the system evolves, the logarithmic negativity becomes zero for any d [r] = —0.9, the right 
panel) and we can say that the system becomes separable at this time. This behavior is not 
changed for the other value of the region size n. The spatial structure of entanglement for 
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this lattice model is simple and we only pay attention to the behavior of entanglement for 
d = 0. 

Then, we consider the evolution of the entanglement. Fig. [3] shows the evolution of 
symplectic eigenvalues z/„ and z/_ for d = with n = 4. During the entire period of 
evolution, the value of is greater than 1/2 and the physical condition fllOp is always 
satisfied. On the other hand, the value of z/_ is smaller than 1/2 initially, then increases and 
exceeds the value 1/2. Thus, the initial entangled state changes into the separable state. 




-2.0 



-1.5 



-1.0 



-0.5 



FIG. 3: Evolution of symplectic eigenvalues for d = 0, n = 4. The upper line represents and 
the lower line represents z>_. The physical condition j/_ > 1/2 is always satisfied. 
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FIG. 4: Ej\f{d = 0) as a function of the conformal time r] for n = 4 case. After rj = rjc ^ —1, Ej\i 
becomes zero and the system is separable. 



We interpret these symplectic eigenvalues behaviors using the logarithmic negativity. Fig. H] 
shows the logarithmic negativity for ci = as a function of conformal time. At some critical 
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time ~ — 1, the logarithmic negativity En becomes zero and the initially entangled state 
changes into a separable state after rjc. As the quantum state, we assume the Bunchi-Davis 
vacuum which imposes the Minkowski vacuum state in the short wave length limit. Thus, 
the entanglement of the scalar field before rj^ implies the remnant of the entanglement of 
the Minkowski vacuum. After 77 = rj^, the regions A and B do not have quantum correlation 
and we expect that the correlation between two regions can be mimicked by an appropriate 
classical distribution function. To understand what time scale determines the critical time 
?7c, we observed how the critical time t]c varies when we change the region size n. 


-5 
-10 
-15 
-20 

FIG. 5: The relation between the critical time rjc and the region size n. The deviation from a linear 
relation represents the effect of boundary condition. 

Fig. shows the relation between the critical time rjc and the region size n. For small value 
of n, the relation coincides with the line rj^. = —n (dotted line). By restoring the dimension 
of the variables, this relation corresponds to 

a{ric) X nAx = —. (36) 
H 

Thus, the quantum entanglement between the regions A and B disappears when the physical 
size of each region exceeds the Hubble horizon length. We must keep in mind that the 
physical size of each region grows as a{r]) x nAx by the cosmic expansion. For larger values 
of n, the relation deviates from the line r]c = —n and we suppose this is due to the effect 
of the periodic boundary condition we imposed for the lattice model. The appearance of 
separability or disentanglement at the Hubble horizon scale coincides with our naive classical 
picture of the separate universe; the Hubble size inflationary domains evolves independently 
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and they can be treated as independent Friedman- Robert son- Walker universes. As is shown 
here, the quantum entanglement between the Hubble size regions is lost and these regions do 
not have the quantum correlation. This is quantum version of the separate universe picture. 
Hence, we can expect that the quantum fluctuations in these regions behave classical. 



IV. ENTANGLEMENT AND QUANTUM TO CLASSICAL TRANSITION 

To clarify the condition of quantum to classical transition of fluctuations in the inflation- 
ary universe, we consider the relation between the disentanglement and the classicalization 
of quantum fluctuations. We can say that the quantum fluctuations become classical if they 
are mimicked by appropriate classical stochastic variables and their statistical nature cannot 
be distinguished from quantum ones. More precisely, the bipartite quantum system under 
the consideration is deflned to be classical if there exists a positive normalizable distribu- 
tion function V on the phase space and the following relation holds for any function F of 
canonical variables: 

{F{qA,PA,qB,PB)) = j dqAdpAdqBdpB'P{qA,PA,qB,PB)F{qA,PA,qB,PB)- (37) 

The left-hand side of this relation is evaluated with respect to the considering quantum state. 
If this relation holds, any correlation function of quantum variables in the bipartite system 
can be imitated by the classical distribution function V on the phase space and appropriate 
classical stochastic variables. 



As is shown by Simon 1 161] and Duan et a/.[17|, for a bipartite two mode Gaussian system, 



if the system is separab 



of P representation: 2l| 



e then the state of the system can be written by the following form 



p = J (fad^p P{a, /3) |a, f3) {a, /3| (38) 

where \a,P) = is the product of the coherent state of A and B, and P{a,P) is a 

positive normalizable function called P function. For a Gaussian state with the covariance 
matrix V, the P function is also a Gaussian function and given by 

p(l) = ±^VMPexp (^-^fP^) ' p=(v + ^ns^sn^^ . (39) 
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where S G Sp(2, R) Sp(2, R) is the symplectic transformation that transforms the covari- 
ance matrix V to the following standard form jl?! : 



Vjj = SVS^ 



/ ar cr \ 

a/r c'/r 

cr ar 
\ c'/r a/r J 

From the definition of P function (l39l). we have 



a — 


d 


J a- 


c 



(40) 



P = 5^ ( Vii 



-1 



(41) 



and the existence of a positive normalizable P function is guaranteed by the condition 
{Vii — -f/2) > 0. In terms of the components of Vn, this condition is represented as 

1 



(a — \c\){a — \c\) > 



(42) 



On the other hand, the symplectic eigenvalues of V, which is invariant under symplectic 
transformations, are 



9^ 



a — 


|c|)(a- 




(a + 


\c\){a + 


\c'\), 


cc' 


> 


a — 


\c\){a + 




(a + 


\c\){a- 




cc' 


< 


a — 


\c\){a + 




(a + 


\c\){a- 




cc' 


> 


a — 


\c\){a- 




(a + 


\c\){a + 




cc' 


< 



(43) 



(44) 



Therefore, the condition of existence of positive normalizable P function is equivalent to the 
separability condition (fTTj) provided that the physical condition (ITOl) is satisfied. From the 
definition of P representation fl55]) . if the system is separable, it is possible to calculate the 
quantum expectation value of the normally ordered product of any operators using the P 
function as a distribution function 



: F{qA,PA,qB,PB) ■) = I dqAdpAdqBdpBP{qA,PA,qB,PB)F(qA,PA,qB,PB) 



(45) 

However, the existence of P function is not sufficient for the establishment of classicality of 
the system; it only guarantees the existence of the distribution function for the normally 
ordered quantities. 
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To derive the condition of classicality of the quantum field, we introduce a Wigner dis- 
tribution function on the phase space and its form for a Gaussian state is given by 

1 /I 



The normahzable Winger function exists for V > and this condition is weaker than the 
physical condition of the state ([5]). Hence, there exists a normalizable Wigner function, 
which does not represent the physical state. The Wigner function ( H6l) for V > is positive 
normalizable and can be interpreted as a distribution function giving the expectation value 
for the symmetrically ordered product of operators 

{{F{qA,PA,qB,PB)}syJ = J dqAdpAdqBdpBW{qA,PA,qB,PB)F{qA,PA,qB,PB)- (47) 

If the difference between the P function and the Wigner function is negligible, these distri- 
bution functions return the same answer for the expectation value of any operator F and 
we have the relation 

{■■F:)^ ^ (F). (48) 

This means that noncommutativity between operators is negligible and the P function or 
the Wigner function plays a role of the classical distribution function, which reproduces the 
quantum expectation value for any operators. Hence, the condition of classicality (1371) is 
established. 

We look for the condition for the establishment of the relation (148|) . For this purpose, 
it is sufficient to consider the condition for the standard form of the covariance matrix Vu 
because this form of the covariance matrix is related to the original covariance matrix V 
via a symplectic transformation. In terms of the covariance matrix, the condition for the 
classicality P ^ W is given by 

Vj^^{Vn-I/2r'. (49) 
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We write down the components of Vjj and {Vjj — 1/2) ^ explicitly 

/ a/r c/r \ 

a^ — c^ 

ar 

c/r 



a^ — c^ 



a^ — 



a^ — c^ 



c'r 



a/r 
a^ — 



c'r 
a^ — c'2 



v 

( {a-l/2r)/r 
(a-l/2r)2-c2 

c/r 



a^ — c'2 



ar 

a2 - c'2 / 



c/r 



(a - l/2r)2 - c2 
(a-r/2)2-c'2 ~(^_^/2)2_c'2 



(a — r/2)r 



c r 



(a-l/2r)2-c2 



c'r 



(a - l/2r)/r 
(a - l/2r)2 -c2 



v 



(a - r/2)2 - c'2 



(a — r/2)r 
(a - r/2)2 - c'2 / 



Thus, if the condition 

<• _ 2 2 2 

ji = a r = a 



a — 


|c'| 


a — 


c 



> 7, /2 = ^ = 



a — 


c 


a — 


|c'| 



1 

»4 



(50) 



4 r^ 

is satisfied, Vjj ~ (V// — //2)^^ and the P function equals the Wigner function. We rewrite 
the condition ( 150|) in terms of the symplectic eigenvalues. For cc' > 0, z/^ = (a — |c|)(a— |c'|) 
and 



/i 



a2iy2 



{a-\c\y 



a — c 



m2 > 



and the condition (|50|) is satisfied if i/_ 1/2. For cc' < 0, 



/i 



a2i/2 



>^!, /2 



a2^;2 



a — c 



./h2 



(51) 



(52) 



and the condition ( l50i) is satisfied if z/_ ^ 1/2. Combining these two cases, we have the 
following result: 



1 1 

^-»2' '^-»2 



V7/ ^ (F,, - //2)-i 



(53) 



Therefore, z/_ ^ 1/2 is the sufficient condition for the system can be treated as classical. 

We can check whether this condition is satisfied in our lattice model. As is show in Fig. [3l 
before the critical time r^^^— 1, z^->l/2 and z/_ < 1/2 and the system is entangled. After 
rjc-, the value of z/_ becomes greater than 1/2 and increases in time. In our lattice model, 
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the relation z/_ < z/_ always holds. For — 1 ^ 77 < 0, the behavior of V- is approximately 
given by 

(^>_)2 _ 1 ~ ^^ll^ _ 0.115144 + 0(7^2) (54) 

and after rj = rjc, the condition z/_ ^ 1/2 is rapidly realized. Hence, the difference between 
the P function and the Wigner function becomes negligible in one Hubble time after the 
system becomes separable at the horizon crossing. As a subset of the separability condition 
for Vji, the inequality 

a' > \ (55) 

holds for the physical state and this corresponds to a standard uncertainty relation. The 
condition of classicality z/_, z>_ ^ 1/2 leads to ^ 1/4 and this also implies the noncom- 
mutativity between canonical variables can be neglected when we evaluating the expectation 
values of operators. This is consistent with the result obtained in the paper l|, 5, 6, 7, 8]; for 
superhorizon scale quantum fluctuations, the noncommutativity between canonical variables 
becomes negligible because the growing mode solution is dominant. In other words, we can 
neglect h in the uncertainty relation. We derived the equivalent condition for the classicality 
from the condition of the existence of the classical distribution function and the symplectic 
eigenvalues. 



V. SUMMARY AND CONCLUSION 



We investigated the appearance of the classical distribution function for the quantum 
fluctuation in the inflationary model using the lattice model of the scalar field. By following 
the evolution of entanglement between two spatially separated regions, we found the clas- 
sicality of the quantum field appears as follows: Initially, when the size of the considering 
region is smaller than the Hubble horizon, the quantum field is in the entangled state. As the 
Universe expands, the quantum state becomes separable when the size of the region equals 
the size of the Hubble horizon. At this stage, the quantum correlation between neighboring 
regions is lost. Then, within about one Hubble time after the horizon crossing, noncom- 
mutativity of operators becomes negligible and the system can be treated as classical. Any 
quantum expectation values can be evaluated using the P function or the Wigner function. 
In other words, there appears classical stochastic nature for variables, which mimics the 
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original quantum dynamics. As we have shown, disentanglement is not a sufficient condi- 
tion for the establishment of classicality of fluctuations. This condition only guarantees the 
loss of EPR-type nonlocal correlations, which are peculiar to quantum mechanics. 

In our analysis, we defined a bipartite system as the subsystem of the entire universe 
and we discard the unobserved dynamical degrees of freedom outside of the observed region. 
Thus, our bipartite system evolves in a nonunitary way and this definition of our system 
effectively takes into account the decoherence mechanism of the considering region. 

We comment on the relation of our analysis to the stochastic approach of inflation 22| 
which treats the quantum dynamics of infiaton fields as the classical stochastic process. 
By coarse graining the scalar field on the large scale {eH)^^,e <^ 1, it can be shown that 
the coarse grained field obeys the Langevin equation and the dynamics of the quantum 
infiaton field is replaced by the classical stochastic process. In the stochastic approach, the 
classical nature of the infiaton field is guaranteed by the appropriate small value of the coarse 
graining parameter e. However, in this approach, the connection between the probability 
distribution and the state of the infiaton field is not clear. The stochastic approach assumes 
the existence of the classical probability distribution from the ffist. However, as we have 
shown in this paper, it is possible to define a probability distribution function only when the 
system becomes separable. We expect that the condition e <^ 1 corresponds to the condition 
z> ^ 1/2 which is stronger than the separability. Anyways, we must reconsider the meaning 
of the probability in the stochastic approach from the view point of entanglement. We will 
report on this topic in a separate publication. 

In this paper, we assumed the Bunch-Davis vacuum state. Previous analysis by J. Les- 
gourgues et a/.|7| considered the nonvacuum initial states that are non-Gaussian and con- 
cluded that the non-Gaussian nature of the state does not affect the transition to the classical 
behavior. However, from the viewpoint of the entanglement, the condition of the separa- 
bility for non-Gaussian states is unknown and the determination of a classicality condition 
for such states is an unsolved problem. Further, we considered the quantum to classical 
transition based on the bipartite entanglement only. This is because the criterion on the 



separability for the general A^-partite system is unknown 1 151]. However it is necessary to 
look for the classicality condition for non-Gaussian states and the A^-partite system to fully 
understand the mechanism of classical to quantum transition of primordial fluctuation. This 
is a future problem to be tackled. 
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